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ABSTRACT 


Helioseismic frequency observations provide an extremely accurate window 
into the solar interior; frequencies from the Michaelson Doppler Imager (MDI) on 
the Solar and Heliospheric Observatory (SOHO) spacecraft enable the adiabatic 
sound speed cand adiabatic index T, to be inferred with an accuracy of a few parts 
in 10 4 , and the density p with an accuracy of a few part s in 10 1 . This has become a 
serious challenge to theoretical models of the Sun. Therefore, we have undertaken 
a self-consist ent,, systematic study of the sources of uncertainties in the standard 
solar models. We found that the largest effect on the interior structure arises from 
the observational uncertainties in the photospheric abundances of the elements 
(e.g., 15% uncertainties in C, N, 0, and Ne, leading to a ~ 10% in the Z / A ratio), 
which affect the sound speed profile at the level of 3 parts in 10 b The estimated 
4% uncertainty in the OPAL opacities could lead to effects of 1 part in 10' 1 ; the 
~ 5% uncertainty in the basic pp nuclear reaction rate would have a similar effect, 
as would uncertainties of ~ 15% in the diffusion constants for the gravitational 
settling of helium. The ~ 50% uncertainties in diffusion constants for the heavier 
elements would have nearly as large an effect. Different observational methods for 
determining the solar radius yield results differing by as much as 7 parts in 10 4 ; 
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we found that this leads to uncertainties of a few parts in 10 3 in the sound 
speed in the solar convective envelope, but has negligible effect on the interior. 
(We did not. explicitly consider the effects of rotational mixing or uncertainties 
in the interior equation of state, which other investigators have found to yield 
uncertainties in the sound speed of order a part in 10 3 .) Our reference standard 
solar model (with Z/X — 0.0245) yielded a convective envelope position R ce = 
0.7135 /?(.,, in excellent agreement with the observed value of 0.713 ± 0.001 R^, 
and was significantly affected (±0.003 R G ) only by Z/X, the pp rate, and the 
uncertainties in helium diffusion constants. Our reference model also yielded 
envelope helium abundance Y e = 0.2424. in good agreement, with the range 0.24 < 
Y, < 0.25 inferred from helioseismic observations; only extreme Z/X values 
yielded Y c outside this range. We found that other current uncertainties, namely, 
in the solar age and luminosity, in nuclear rates other than the pp reaction, in 
the low-temperature molecular opacities, and in the low-density equation of state, 
have no significant effect on the quantities that can be inferred from helioseismic 
observations. The predicted pre-main-sequence lithium depletion is uncertain 
by a factor of 2. The predicted neutrino capture rate is uncertain by ~ 30% 
for the 37 C1 experiment and by ~ 3% for the 71 Ga experiments (not including 
uncertainties in the capture cross sections), while the 8 B neutrino flux is uncertain 
by ~ 30%. 

Subject headings: diffusion — neutrinos — Sun: abundances — Sun: helioseis- 
mology — Sun: interior 


1. Introduction 

Helioseismology provides a powerful tool to explore the deep interior of the Sun. Mea- 
surements of solar interior quantities are provided with unprecedented accuracy. The Michel- 
son Doppler Imager (MDI) on the Solar and Heliospheric Observatory (SOHO) spacecraft 
(Rhodes et al. 1997), the Global Oscillation Network Group (GONG), the Birmingham So- 
lar Oscillation Network (BiSON: Chaplin et al. 1996), and the Low -l instrument (LOWL: 
Tomczvk et. al. 1995) provide helioseismic frequency measurements with accuracies of a few 
parts in 10 5 . From these, the sound speed c throughout most of the solar interior can be 
inferred with an accuracy of a few parts in 10 4 , as can the adiabatic index Tp the density 
can be inferred with an accuracy of a few parts in 10 3 (Basu, Pinsonneault, & Bahcall 2000; 
Bahcall, Pinsonneault, & Basu 2001). 

The ultimate goal of our work was to explore systematically a wide range of solar models 
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with relatively modest, amounts of mass loss on tin* early main sequence, and to test theii 
viability via helioseisrnological measurements (as well other observational constraints). This 
mass loss investigation is presented in our companion paper "Our Sun \ (Sackmann <!c 
Boothroyd 2001). Since the consequences of moderate early solar mass loss are expected to 
be small, it is essential to understand the consequences of uncertainties in the input physics 
and in the input parameters of the solar model. The physics inputs include the equation of 
state, opacities, nuclear rates, diffusion constants, the treatment of convection (including the 
possibility of overshoot), and the effects of rotation and mass loss. Observed solar parameters 
include the solar age and the present solar radius, luminosity, surface composition, and solar 
wind. A considerable amount of work has been published investigating many of the above 
effects, as will be discussed in § 3. However, some basic uncertainties still warrant further 
attention. Before proceeding to our mass loss work, we found it necessary to try to extend the 
investigations of the above uncertainties; in particular, the major consequences arising from 
uncertainties in the present observed solar surface Z/X ratio and in the basic p-p chain rate 
have not been sufficiently considered in the recent (most accurate) work of other investigators. 
To determine the consequences of these uncertainties, we computed st andard solar models 
with various values of the Z/X ratio and the p-p chain rate (lying in the permitted ranges). 
To obtain a set of self-consistent results, we also investigated the consequences of a number 
of other uncertainties, namely, in other nuclear rates, in envelope opacities, in the equation 
of state, in diffusion constants, from different methods of handling interior opacities, and in 
the solar age. It is the aim of this paper to present the effects of the above uncertainties on 
the run of the sound speed and density in the solar interior, and on the radius R cc of solar 
convection and the solar envelope helium abundance I, We also present, the effects on the 
solar lithium depletion (excluding rotation effects) and on the production of solar neutrinos. 
(Note that lithium depletion in a non-rotating standard solar model occurs entirely on the 
pre- main-sequence; this predicted depletion is significantly smaller than the observed solar 
lithium depletion factor. Rotational mixing on the main sequence is generally invoked to 
explain the remaining lithium deficit). 


2. Methods 

We computed a reference standard solar model using up-to-date physics and observed 
solar parameter values, and several dozen variant solar models in which one of these “inputs 
was varied within the allowed errors. (In a few cases, more than one of the “inputs” was 
varied, or the size of the variation exceeded the size of the quoted errors in order to get a 
better estimate of the sensitivity.) By comparing with the reference standard solar model, 
the sensitivity to the uncertainties in the “inputs” could be determined. The stellar evolution 
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code used to compute these solar models is descended from that used earlier in our “Sun 111" 
paper (Sackmann, Boothroyd, k Kraemer 1993; see also Boothroyd k Sackrnann 1999), 
but has been extensively updated for improved accuracy, including provision for much finer 
zoning as well as up-to-date input physics. 

The reference standard solar model used the OPAL equation of state (Rogers, Swenson, 
k Iglesias 1996) in the interior, and the MHD equation of state (Dappen et al. 1988) where it 
was designed to be accurate, namely, in the outer envelope at log p < —2 (this corresponds 
to r > 0.94 Rq and log T < 5.5 in the present Sun). A version of the MHD equation of state 
computation program was kindly provided to us (D. Mihalas 1999, private communication), 
and minor modifications allowed computation of the MHD equation of state for various hy- 
drogen abundances and metallicities down to pre-main-sequence photospheric temperatures 
and even below (A. I. Boothroyd 2001, in preparation). In both cases, the equation of state 
was interpolated in metallicity as well as in hydrogen abundance, temperature, and density, 
in order to take into account metallicity variations due to diffusion and nuclear burning. 
Variant cases were tested where the OPAL equation of state was used down to log T = 4 or 
even all the way down to its lower limit of validity at log T — 3.7 — in the lat ter case, the 
MHD equation of state was relevant only to the pre-main-sequence evolution. The MHD 
equation of state tables computed for this work considered H and He, plus a 13-element sub- 
set of the Grevesse k Noels (1993) heavy element composition (i.e., C, N, 0, Ne, Na, Mg, 
Al, Si, S, Ar, K, Ca, and Fe); besides the effects due to neutrals and ions, effects due to H _ , 
H 2 , and HJ were accurately taken into account, and approximate effects of the molecules 
C 2 , N 2 , 0 2 , CH, CN, CO, NH, NO, OH, C0 2 , and H 2 0 (although these molecules have no 
effect for solar models: they were actually added in anticipation of use in asymptotic giant 
branch models). The OPAL equation of state considers only H, He, C, N, O, and Ne (plus 
hydrogen-molecule effects); however, the equation of state is quite insensitive to the precise 
makeup of the metallicity, and the OPAL equation of state may well be the more accurate 
one in the region where both are valid (Rogers et al. 1996). 

The reference standard solar model used the 1995 OPAL opacities “s:opaL:GN<m” (Igle- 
sias k Rogers 1996) at temperatures log T > 4, and the Alexander k Ferguson (1994) 
opacities “/tAiexander” (which include molecular opacities) at log T < 4. Variant cases tested 
the use of the Sharp (1992) molecular opacities “/csharp” instead at log T < 4, and other 
high-temperature opacity cases were also considered. The standard 1995 OPAL opacities 
use the Grevesse k Noels (1993) heavy element composition for the components of the 
metallicity Z, but the online opacity computation feature of the OPAL web page 2 allowed 
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computation of OPAL opacities appropriate to the more recent Grevesse Sauval (199S) 
mixture (“kopaL:GS 89”); these were tested in variant models, as were opacities appropriate 
to the older Grevesse (1984) mixture (“kopaliGiW’)- The even older and much less accurate 
Los Alamos (LAOL) opacities “/claolss” (Keadv 1985, private communication) were also 
tested in one case. 

In our reference standard solar model, we did our best to account for temporal and 
spatial variations in the opacity due to composition changes from diffusion and nuclear 
burning. The OPAL opacity tables allow interpolation of the opacity as a function of the 
hydrogen abundance and the metallicity Z; the abundances of the metals comprising Z 
are always proportional to Z in the OPAL tables, i.e., a “scaled solar distribution. In 
addition, the OPAL opacity tables contain mixtures with excess carbon and oxygen (beyond 
that contained in the scaled solar metallicity), allowing interpolation in carbon and oxygen 
abundances. 

To take into account the variations in abundance due to diffusion, the metallicity 
value Z K that we used for metallicity interpolation in the OPAL opacity tables was pro- 
portional to the abundances of the elements heavier than oxygen: Z K = Z/,, where Z h = 
z o[Hheavy X i}/{T / heavy( X i)o]' where Z o and Mo are the protosolar metallicity and com- 
position, respectively, and J2 heavy refers to a sum over elements heavier than oxygen. In 
other words, we scaled the initial solar metallicity by the shift in heavy element abundances 
resulting from diffusion (note that our diffusion routines assumed that all metals diffused 
alike). However, for the CNO elements, there are additional changes due to nuclear burning, 
so that the CNO abundance profiles are not proportional to the heavy element abundance 
profiles; Turcotte et al. (1998) find that conversion of C and O into N results in an opac- 
ity change of ~ 1% that cannot be modelled by a change in Z alone. These variations in 
the CNO elements relative to Z/, were accounted for in an approximate manner by an addi- 
tional two-dimensional interpolation in nominal “excess carbon and oxygen” abundances C PX 
and 0«, where these excess abundances account for any variation in the CNO abundances 
relative to the scaled solar metallicity Z K of the OPAL opacity tables. Since there were 
no explicit opacity tables for changes in nitrogen, the best one could do was to distribute 
excess nitrogen equally between C eT and O ex (i.e., to assume that nitrogen opacities were 
midway between those of carbon and oxygen): C ei = C — C 0 Z K /Z 0 -L 0.5(N — N 0 Z K /Z 0 ) and 
O ex = O - O 0 Z K /Z 0 + 0.5(N - N 0 Z*/Z 0 ), such that CO ex = C ei + O cx = Z - Z K . 

Note that in general, as first carbon and then oxygen is burned to nitrogen, either C ei 
or O ex will be negative, implying extrapolation of the OPAL tables in the direction of zero 
C or O by a non-negligible fraction of Z. One might consider the uncertainties in such an 
extrapolation to be worse than the error inherent in treating all CNO opacities alike; in this 
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case, if one of C ex or O ex was negative, one would set it to zero and subtract an equivalent 
amount from the other — this was the method followed in our reference models. However, 
we also tested the case where negative values were allowed, extrapolating as required. These 
turned out to yield essentially identical results to the reference model: the rms difference 
was only 0.00006 in the sound speed profile A cjc and 0.0003 in A p/p, comparable to the 
estimated numerical accuracy of the models. An alternative case where this latter CNO 
variation was not approximated by excess CO at all, i.e., having Z K = Z h but CO fI = 0.0. 
also yielded almost identical results (rms difference of 0.00005 in A cjc and 0.0004 in A p/p). 
Similarly negligible rms differences were found for several cases testing different prescriptions 
for defining Zh, C ex , and O ei (differing in whether mass fractions or number densities were 
used, in how the excess nitrogen was divided between C ex and O ex , and in whether negative 
values were allowed for C ei and O ei ). 

An estimate for the upper limit of the effects of opacity uncertainties was made by 
simply setting Z K = f Z for a constant factor / = 0.9, 0.95, 1.0, 1.05, or 1.10; note that the 
case / = 1.0 is only slightly different from the reference standard solar model. 

The reference standard solar model used the NACRE nuclear reaction rates (Angulo 
et al. 1999), supplemented by the 7 Be electron capture rates of Gruzinov & Bahcall (1997) 
for log T > 6 and of Bahcall & Moeller (1969) for log T < 6 (note that this latter low- 
temperature region is irrelevant for solar models). Variant models tested cases with nuclear 
rates changed according to the upper or lower limits quoted in the NACRE compilation; 
one case tested the use of the older Caughlan & Fowler (1988) nuclear rates. The pro- 
gram uses the minimum of the weak (Salpeter 1955), intermediate (Graboske et al. 1973), or 
strong (Itoh et al. 1979; Ichimaru & Utsumi 1983) screening factors; for solar conditions, this 
means weak screening, which is a very good approximation to the exact quantum mechan- 
ical solution (see, e.g., Bahcall, Chen, & Kamionkowski 1998b; Gruzinov & Bahcall 1998). 
Deuterium was not considered separately (it was assumed to have been entirely burned to 
J He on the early pre-main-sequence), but all the other 15 stable isotopes up to and including 
,8 0 were considered in detail (i.e., nuclear equilibrium was not assumed for any of them). 
The other stable isotopes up to 28 Si were included in the code (plus a few long-lived unstable 
isotopes, as well as Fe and a category for the sum of the other elements heavier than Si), 
but their nuclear reactions were not included since there are no significant effects under solar 
conditions (except for 19 F, which was assumed to be in CNO-cycle nuclear equilibrium for 
nuclear rate purposes). Neutrino capture cross sections were taken from Bahcall & Ulrich 
(1988), except for the 8 B-neutrino cross section for capture on 37 C1, where the more recent 
value (5% higher) of Aufderheide et al. (1994) was used. 



A set of subroutines 3 were kindly provided to us (M. H. Pinsonneault 1999, private 
communication) that take into account the diffusion (gravitat ional settling) of helium and 
heavy elements relative to hydrogen (see also Thoul, Bahcall, Loeb 1994, Bahcall, Pinson- 
neault, & Wasserburg 1995). These subroutines assume that all heavy elements diffuse at 
the same rate as fully-ionized Fe; this yields surprisingly accurate results, as may be seen bj 
comparing to the results of a more detailed treatment (Turcotte & Christensen-Dalsgaard 
1998; Turcotte et al. 1998). The upper limit of the effects of uncertainties in the diffusion 
constants was estimated by simply multiplying the diffusion constants, either for helium or 
for the heavy elements, by a constant factor. 

The Schwarzschild criterion was used to define convective boundaries; no core over- 
shooting or envelope undershooting was allowed. Note that Morel, Provost, X Berthomieu 
(1997) found that including convective core overshooting had a negligible effect on the solar 
sound speed and density profiles, but on the other hand that including convective envelope 
undershooting bv even a tenth of a pressure scale height moved the solar convective envelope 
boundary inwards by eight times the uncertainty in the observed value, yielding a sharp spike 
in the difference between observed and calculated sound speed profiles. Rotation-induced 
mixing was not considered in our models; the effect that it would have is discussed in § 3.1. 

The reference standard solar model used a value of Z / A — 0.0245 for the the piesent so- 
lar surface metals-to-hydrogen ratio (by mass fraction), as given by Grevesse & Noels (1993). 
Variant models tested a 10% lower value of ZjX — 0.022, and the ~ 13% higher older value 
of ZjX - 0.0277 (Grevesse 1984; Anders & Grevesse 1989). Additionally, variant models 
with the value ZjX — 0.023 recommended by Grevesse & Sauval (1998) were considered. 
The presolar abundances of the heavy elements were taken to be the ones used in the OP AL 
opacities used in the solar model; for the reference standard solar model and most variants, 
this was the mix of Grevesse & Noels (1993). A variant model with ZjX - 0.023 also tested 
use of the heavy element mix and OPAL opacities appropriate to the abundances reported 
by Grevesse & Sauval (1998), and a variant model witli ZjX = 0.0277 tested use of the 
heavy element mix and OPAL opacities appropriate to the abundances reported by Grevesse 
(1984). Tests were made with C, N, O, and Ne abundances increased or decreased by their 
uncertainties of 15% relative to Fe (with the OPAL opacities appropriate to these revised 


A present solar mass of M & = 1.9891 x 10 33 g (Cohen &’ Taylor 1986) was used in all 
cases. The uncertainty of 0.02% is too small to have any significant effect, and is in fact 
smaller than the amount of mass lost by the sun since it was formed in the form of radiation 
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alone (photons and neutrinos: AM = AE/c 2 ). Mass loss from this effect and from the solar 
wind was ignored for all cases considered in this paper; as is shown in our companion paper 
“Our Sun V” (Sackmann k Boothroyd 2001), it would have negligible effect, even if the 
average solar wind over the past 4.6 Gyr had been an order of magnitude higher than its 
present value, as is suggested by measurements of noble gas isotopes in lunar rocks (Geiss 
1973; Geiss k Bochsler 1991; Kerridge et al. 1991). 

For the reference standard solar model and most variant cases, a present solar luminosity 
of Lq = 3.854 x 10 33 erg s _1 was used, as discussed in Sackmann et al. (1993). A value 
0.3% lower (3.842 x 10 33 erg s -1 ) with an estimated 1 -a uncertainty of 0.4% was recently 
obtained by Bahcall et al. (2001), based on the observations of Frohlich k Lean (1998) and 
Crommelynck et al. (1996). Variant models considered the effect of using this more recent 
solar luminosity value, and high and low values 2-cr (0.8%) above and below it. Note that 
most of the uncertainty in L 0 comes not from uncertainties in the present solar irradiance, 
but rather from uncertainties in the slight long-term variability of the solar luminosity. For 
example, Lean (2000) estimates a difference of 0.2% between the present value and that 
during the seventeenth century Maunder Minimum, and even this timescale is much shorter 
than the Sun’s Kelvin-Helmholtz timescale of ~ 10 Myr. 

The reference standard solar model used a solar radius at the photosphere (r = 2/3) of 
Rq — 695.98 Mm (Ulrich & Rhodes 1983; Guenther et al. 1992). Variant models considered 
the effect of using the value of 695.78 Mm suggested by the helioseismic /-mode study of 
Antia (1998), or the value of 695.508 Mm suggested by the solar-meridian transit study of 
Brown k Christensen-Dalsgaard (1998). 

In his Appendix to Bahcall et al. (1995), G. J. Wasserburg provides a systematic analysis 
of the upper and lower bounds on the age of the Sun, as obtained from isotopic ratios 
measured in meteorites. We paraphrase his discussion in this paragraph: The protosolar 
nebula (out of which the Sun and the meteorites formed) contained not only the stable 
isotope 27 Al but also the unstable isotope 26 Al, which decays into 26 Mg with a half-life of 
only 0.7 Myr. This 26 A1 must have been injected into the protosolar nebula from the stellar 
source where it was created. Isotopic measurements of meteoritic crystallized refractory 
condensates show that they had a ratio 26 Al / 27 Al = 5 x 10~ 5 at the time they formed. Even 
if the stellar source of the 26 Al had a very high ratio 26 A1/ 27 A1 ~ 1. the time interval between 
the formation of the 26 Al and the formation of the meteorite cannot have exceeded ~ 11 Myr. 
The Sun cannot have formed later than these meteorites. However, the Sun/meteorite system 
must have formed after the injection of 26 Al into the protosolar nebula, so the Sun cannot 
have formed earlier than ~ 11 Myr before the formation of these meteorites. The age of the 
meteorites has been accurately measured using 207 Pb/ 206 Pb ratios, to be 4.565 ± 0.005 Gyr. 



9 - 


It follows that the Sun cannot have formed more than 4.59 Gyr ago, nor later t han 4.o5 Gyi 
ago (i.e., 4.565 + 3 x 0.005 -t- 0.011 Gyr. or 4.565 — 3 x 0.00a Gyi). 

Our solar models were started relatively high on the pre-main-sequence Hayashi track, 
with central temperatures below 10 c K; note that the ages t e of all our models are quoted 
relative to this pre-main-sequence starting point, which is within a few Mvr of the solar 
formation age constrained by the meteoritic ages discussed above. A relatively high value 
of f 0 = 4.6 Gyr was used for the reference standard solar model; to get a reliable estimate 
of the sensitivity of the models to the solar age uncertainty, variant models were computed 
with ages differing by very large amounts, namely, t© = 4.5 and 4.7 Gyr. 

From our Hayashi track starting point, it takes only ~ 3 Myr for the luminosity on the 
pre-main-sequence to drop below 1 L© (i.e., below the present solar luminosity), but much 
longer, namely an additional 40 Myr, to reach the zero-age main sequence (ZAMS) — we 
have defined the ZAMS as the stage where the pre-main-sequence contraction terminates and 
the Sun begins to expand slowly, as nuclear burning in the core (rather than gravitational 
contraction) supplies essentially all of the solar luminosity. For the next ~ 50 Myr, evolution 
is fairly fast, as p + C reactions burn the initial carbon to nitrogen near the Sun’s center, 
resulting in a short-lived convective core. Subsequently, the carbon abundance approaches 
its CN-cycle equilibrium abundance, the convective core disappears, and the Sun settles 
down to burn hydrogen mainly via the pp- chain reactions. Note that the pre-main-sequence 
timescale implies that the total solar age t & used in this paper can be converted into a main 
sequence solar lifetime by subtracting about 0.04 Gyr. 

We compared our solar models to profiles of the solar sound speed c©, density p©, and 
adiabatic index (r,)© obtained from the helioseismic reference model of Basu et al. (2000) 4 . 
which they obtained by inversion from the helioseismic frequency observations. In the inver- 
sion process, a standard solar model is required, but Basu et al. (2000) demonstrated that 
the resulting c© and p© profiles of the helioseismic reference model are relatively insensitive 
to uncertainties in the standard solar model used for this purpose (except for uncertainties 
in R e , as discussed in § 3.1). They estimated a net uncertainty of few parts in 10 4 for 
the sound speed c© and adiabatic index (Ti)©, and a few parts in 10 for the density p©. 
However, in the Sun’s core (r < 0.1 /£©), systematic uncertainties in the helioseismic sound 
profile are increased by a factor of ~ 5; this was demonstrated by Bahcall et al. (2001), 
who compared helioseismic inversions of different helioseismic data sets. We used their com- 
parison to estimate the r-dependence of the systematic error in c© in the core and in the 
convective envelope (namely, a fractional systematic error decreasing linearly from 0.0013 at 


4 From the denser-grid machine-readable form of their Table 2, at http://www.sns.ias.edu/~jnb/ 
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r - 0.05 R 0 to 0.0003 at r = 0.2 Rq, constant, from there to r = 0.72; 7?©, then increasing 
linearly to 0.00052 at r = 0.94; Rq). For c©, this systematic error can be significantly larger 
than the statistical errors, and we combined the two in quadrature to get the fractional 
error ( o c jc ) for the purpose of calculating weighted rms differences — the rms fractional dif- 
ference in c is given by ({^3 [(Ac/c)/(<r c /c)] 2 } / {13 (V^c/c)] 2 }) ■ For (Ti)© and p Q , the 

systematic errors are comparable to or smaller than the statistical ones, and the statistical 
errors sufficed for calculating weighted rms differences. 


3. Results and Discussion 
3.1. Sound Speed and Density Profiles 

We present in Figures 1 through 5 profiles of the adiabatic sound speed differences 
5c j c = ( Cq - c mode i) / Cq] profiles of the density differences 5p/p = (p© - p mo del)/po > are 
available online 5 . For our equation of state comparison, we also considered the equivalent 
fractional difference in the adiabatic index Fj. The theoretical sound speeds c modet and 
densities p mo dei are from our computed reference standard solar model and from our variant 
standard solar models. Our reference standard solar model used current input parameters, 
as discussed in § 2; our variant standard solar models comprised standard solar models 
with one or more input parameters varied within the permitted range. We present all our 
sound speed and density profiles in terms of differences relative to the observed helioseismic 
reference profiles. This choice of presentation not only allows one to see the effects of the 
uncertainties in the input parameters, but also shows which choice of input parameters agrees 
best with the helioseismic observations. 


Rotation effects: The prominent peak in 5cjc at r ~ 0.7 Rq visible in Figure la is due to 

the neglect of rotation-induced mixing just below the base of the solar convective envelope, 
as has been shown by a number of investigators who have included parameterized rotational 
mixing (see, e.g., Richard et al. 1996; Brun, Turck-Chieze, & Zahn 1999; Basu et al. 2000; 
Bahcall et al. 2001). There are significant uncertainties in the physical processes that lead 
to rotation-induced mixing. However, all of these investigators agree that rotational mixing 
is capable of smoothing out the peak at r ~ 0.7 Rq, and that this has a relatively small 
effect on the sound speed elsewhere in the Sun (a fractional change of less than 0.001). For 
example, Bahcall et al. (2001) found that including “maximal” rotational mixing spread 


5 http://www.krl.caltech.edu/~aib/papdat.html 
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out this peak over the region 0.3 11q < r < 0.7 Rq . eliminating the j)roininent peak but 
worsening the agreement with the helioseismological sound speed profile by about 0.001 in 
much of the solar interior (0.3 Rq < r < 0.6 Rq). The “minimar mixing model of Richard 
et al. (1996) yielded much the same result, as did a similar model of Brun et al. (1999). As 
far as the core and the convective envelope are concerned, rotational mixing should have no 
significantly effect, as shown by the above authors. 

Since the prominent peak at r ~ 0.7 R @ results from the neglect of rotational mixing, 
we did not require agreement in this region between profiles from our theoretical models and 
profiles inferred from the helioseismic observations. Nor did we require agreement in core 
region, since the present helioseismic observations still result in large uncertainties in the 
inferred profiles there; for example, as shown by Bahcall et al. (2001), the use of a diffeient 
helioseismic dataset could remove the the sharp upturn in 6c/c for r < 0.1 Rq in Figure la 
(and even convert it into a downward trend). On the other hand, we aimed for agreement 
in the regions 0.1 Rq < r < 0.6 R G and 0.72 Rq < r < 0.94 Rq, where disagreements are 
due to imperfections in the input physics or uncertainties in the observed solar parameters. 
This is demonstrated by our variant models and by the variant models of Morel et al. (1997) 
and Basu et al. (2000). 

Convergence accuracy effects: The accuracy with which the model is converged to the 

solar radius, luminosity, and ZjX values can affect the sound speed profiles. The extent 
to which this occurs depends on the accuracy of the convergence, and on the sensitivity of 
the profiles to variations in Rq, Lq, and ZjX (these are discussed in detail below). Oui 
convergence accuracy resulted in effects no larger than a few parts in 10 on the sound speed 
in the solar interior (r < 0.6 Rq). In the convective envelope, where the sound speed is quite 
sensitive to Rq , the effect was typically less than a part in 10 4 (this included the reference 
standard solar model), but could be as high as a few parts in 10 4 for the few models with 
the worst convergence in Rq. This was determined not only by the sensitivity considerations 
discussed below, but also by considering models with differing degrees of convergence on Rq, 
Lq, and Z/X. 

Note that the sound speed profile in most of the solar convective envelope (0.72 R & < 
r < 0.94 Rq) is sensitive mainly to the equation of state and to the solar radius, with other 
uncertainties having only a minor effect there, as can be seen by considering this subregion 
in Figures 1 through 5. 

Zoning effects: We investigated the effects of using two different zonings. Our coarse- 

zoned models had about 2000 spatial zones in the model, and about 200 time steps in the 
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evolution from the zero-age main sequence to the present solar age (plus about 800 time 
steps on the pre-main-sequence). Typically, these models were converged to match the solar 
luminosity and radius to about a part in 10 5 , and the solar surface ZjX to a part in 10 4 ; a 
few cases where convergence was slow were nearly 10 times worse. Our fine-zoned models had 
10000 spatial zones and took 1500 main-sequence time steps (plus 6000 pre-main-sequence 
time steps) — a factor of 5 increase in both spatial and temporal precision — and were 
typically converged to better than a part in 10 5 for R e and and a lew parts in 10’ 
for Zf X. We also tested some very coarse-zoned models, with 1000 spatial zones. 100 main- 
sequence time steps (plus 600 pre-main-sequence time steps), and convergence to the solar 
parameters of a few parts in 10 4 . (A coarse-zoned converged solar model took a few hours 
of CPU-time on a high-performance ES40 computer, as compared to about four days of 
CPU-time for a fine-zoned converged model; these times were roughly tripled on a 450 Mhz 
Pentium III PC.) 

Figure la shows that the fine zoning made only a very modest improvement relative to 
the coarse-zoned case, less even than the statistical errors in the sound speed and density 
profiles obtained from hclioseismic inversions. Even the very coai so- zoned test case did not 
do too badly: in the solar interior, it differs from the fine-zoned case by no more than 
0.0004 in the sound speed profile and 0.003 in the density profile (amounts comparable to 
the systematic uncertainties in the helioseismic inversion) the nns differences (over the 
entire Sun) are even smaller, at 0.0003 and 0.002, respectively. The coarse-zoned model 
did about 3 times better still, with rms differences relative to the fine-zoned case of 0.0001 
and 0.0008 for 6c/c and 6p/p, respectively. Zoning changes had no effect on the adiabatic 
index Tj (the coarse-zoned T] curves were not plotted in Fig. lb, since they would be precisely 
superimposed on the fine-zoned curves). Additional tests demonstrated that differ ent. zoning 
always led to the same negligibly small systematic shift in solar interior 6c/ c and 6p/p 
values (although inaccuracies in matching the observed solar surface parameters could lead 
to slightly larger random variations in the convective envelope region r > 0.7 f? 0 ). We 
therefore felt justified in running most of the models with our coarse zoning. Note that Morel 
et al. (1997), with about 1000 spatial zones, 60 main-sequence time steps, and convergence 
to present solar surface parameters of a part in 10 4 (similar to our very-coarse-zoned case), 
claimed a numerical internal accuracy of 0.0005 in the sound speed, similar to what we found. 

Equation-of-state effects: Morel et al. (1997) compared two relatively-current equations 

of state, and found a small but non-negligible effect. They compared the OPAL equation 
of state (Rogers et al. 1996) with the CEFF equation of state (Christensen-Dalsgaard & 
Dappen 1992), finding differences of slightly over 0.001 in the sound speed profile and up 
to 0.01 in the density. (Note that they had set the value of Z eos used in their equation 
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of state to a fixed value of 0.019; however, since the equation of state is only very weakly 
sensitive to Z, this should have only a minor effect on their models, and should not affect 
their comparison of the two equations of state.) 


Guzik & Swenson (1997) presented more extensive comparisons, relative to the OPAL 
equation of state, of both their own SIREFF equation of state and the MHD equation of 
state of Dappen et al. (1988) (they too used a fixed Z eos , of 0.02, but again this should 
not affect the comparison). They likewise found an effect Ac/c < 0.001 at / < 0.95 P. 
(with differences up to 0.004 near the surface, probably due to the different methods used 
to determine ionization effects, which are very large there due to the helium and hydrogen 
ionization zones). They also presented differences between the values of the pressure P, 
specific heat at constant pressure C p , internal energy U , and adiabatic index R. As one 
might expect from the sound speed, differences between equations of state near the solar 
surface (r > 0.95 P©) were relatively large, of order 1%. However, for r < 0.95 P© Guzik 
& Swenson (1997) reported OPAL - MHD differences A P/P < 0.0015, A U/U < 0.002, 
Ar,/r, < 0.0007, and AC, ,/C,, < 0.015 (with AC, ,/C,, < 0.006 for r < 0.9 P,„); the 
OPAL - SIREFF differences were slightly smaller for C p , slightly larger for P and U, and 
much larger (a factor of ~ 3) for R. For r > 0.3 P©, the SIREFF value of R has several 
relatively large “wiggles” (fractional variations ~ ±0.002) relative to either OPAL or MHD; 
but most of the difference in the core, and perhaps some difference in the average trend 
further out, may be due to the fact that SIREFF includes relativistic electron effects, while 
MHD and OPAL do not. Elliot & Kosovichev (1998) estimated that inclusion of relativistic 
effects would reduce the MHD or OPAL value of I ] by a fraction 0.002 at r ~ 0.1 P©, this 
correction growing less with increasing r, to reach 0.001 at r ~ 0.3 P© and zero neai the 
solar surface. They pointed out that such a shift in R for models using the OPAL or MHD 
equation of state would significantly improve the agreement in the solar interior with the 
inferred helioseismic R profile. Certainly, if such a correction were applied to our R curve in 
Figure lb, the model profile would agree with the helioseismic profile within the statistical 
errors for all r < 0.6 P© (recall that decreasing a model quantity shifts the curve upwards 
in the figure). 


The adiabatic sound speed is defined as c = (r 1 P/p) 1/ ' 2 ; changes in the solar ratio of 
P/ p would result not only from changes in the equation of state but also from readjustments 
of the solar structure in response to these changes, so it is not obvious a priori what effect 
relativistic corrections would have on the sound speed. Consideration of the effect from T] 
alone suggest that relativistic corrections might reduce the slope at r < 0.5 P© in the 6c jc 
curve of Figure la. The sound speed differences presented by Guzik & Swenson (1997) for 
their OPAL - SIREFF comparison suggest that this would in fact be the case, and that 
a fractional decrease of order 0.001 in the sound speed Crnodei uear the Sun s center (i.e.. 
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an increase of 0.001 in Sc/c there) would result from relativistic: corrections to the OPAL 
equation of state. 

Richard et al. (1998) compared the OPAL value of r, with the MHD value in the 
convective envelope (0.72 Rq<t < 0.98 R 0 ) resulting from the equation of state with that 
from the MHD equation of state (Dappen et al. 1988), finding that the OPAL equation of 
state appeared to perform slightly better. This is not very surprising; what is perhaps more 
surprising is how well the MHD equation of state does in the solai inteiioi, since it was not 
originally designed to be accurate for p > 10" 2 g cm -3 (i.e., r < 0.94 and log T < 5.5 
in the Sun). We investigated the effect of changing the equation of state only in this outer 
region where both are expected to be valid. 

Our reference standard solar model used the OPAL equation of state in the interior 
regions, switching over to the MHD equation of state in the outer envelope; this switchover 
was performed gradually over the region — 1.5 > log p > — 2 (corresponding to 5.8 > log T > 
5.5 in the present Sun). We compared this reference standard model with two cases where 
the switchover occurred even further out in the envelope: a case OPALeos-midT wheie the 
switchover occurred for 4.0 > log T > 3.9, and a case “OPALeos-lowT where the switchover 
occurred for 3.75 > log T > 3.7. In this latter “OPALeos-lowT” case, the MHD equation of 
state is used only outside the Sun’s photosphere, and thus has negligible effect on the main 
sequence evolution. On the other hand, any artifacts induced by the switchover might be 
smaller in the former “OPALeos-midT” case, since differences between the two equations of 
state are significantly less near log T = 4 than near log T = 3.7. Note that, while the OPAL 
opacity tables are unreliable for log T < 4 due to their neglect of molecular opacities, the 
OPAL equation of state tables include molecular hydrogen effects, and should be accurate 
down to their lower tabulation limit of log T = 3.699 (Rogers et al. 1996). 

Figure 1 shows that, as one would expect, changing the equation of state in the outer 
envelope alone has no effect on the interior, and only a minor effect in that part of the convec- 
tive envelope where the equation of state remains unchanged. (Since the effects were so small, 
we computed fine-zoned cases for this equation-of-state test; these are the ones presented in 
Fig. 1.) The OPAL equation of state appears to yield slightly better Fj values than the MHD 
equation of state in the envelope at r > 0.9 Rq (log p < —1.5). We agree with Guzik &; 
Swenson (1997) in finding maximum differences between these two equations of state in the 
outer solar envelope (at log T ~ 5) of about 1% in pressure P and in x P = (9 In Pjd\n p}ti 
and of 2% in the specific heat at constant volume C v and in xt = (d In P/d In T) p . 

Solar radius effects: Basu et al. (2000) demonstrated that using a solar radius differ- 

ent from the standard value of Rq = 695.98 Mm could have a small but not completely 
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insignificant effect on both the sound speed profile inferred from helioseismic inversions and 
that computed in solar models. They found that using the 0.03% smaller solar radius value 
R q - 695.78 Mm (case “R78”) suggested by the /-mode study of Antia (1998) would re- 
duce the inferred helioseismic sound speed profile throughout the Sun by about the same 
small fraction, namely 0.0003 (nearly independent of position in the Sun); similarly, using the 
0.07% smaller value Rq = 695.508 Mm (case “R508”) suggested by the solar-meridian transit 
study of Brown & Christensen-Dalsgaard (1998) would reduce the inferred helioseismic sound 
speed profile bv 0.0007. On plots such as ours of fractional differences 8c/c = (c & — c miH i r .i)/ c c-'.y 
such a reduction in “c Q " would shift all the curves downwards by the given amounts. This 
shift has not been performed in our diagrams — we only display our c model values relative 
to the cq values of Basu et al. (2000) — but this effect has been included in our quoted 
rms values relative to the Sun. Basu et al. (2000) also calculated that a change in the solar 
radius would result in a non-uniform shift in the inferred helioseismic p Q and (Tj ) 0 profiles, 
by amounts comparable to the statistical errors in these quantities; although these shifts are 
barely significant statistically, in contrast to the shift in c 0 , for completeness their effects 
have been applied to our rms values calculated relative to the helioseismic profiles foi the 
“R78” and “R508” cases. 


Our “R78” and “R508” models show that changing the value of R G has no effect on F, , 
and only in the convective envelope does it affect c and p significantly, with the largest effect 
being near the solar surface. The effect on c is shown in Figure la. For clarity, these R78 
and “R508” curves are omitted at r < 0.72 in the “peak” region, they differ from the 
reference standard solar model by a few parts in 10 4 , and this difference drops to about a 
part in 10 5 for r < 0.6 Rq. In the convective envelope, the “R78” case is an improvement on 
the reference model (reducing the rms error in 6c/c from 0.0007 to 0.0004 when one includes 
the effect of the shift in the inferred helioseismic profiles), but the “R508” case is worse 
(rms of 0.0011). For the “OPALeos-midT” case, reducing the solar radius always worsens 
agreement in the convective envelope (an envelope rms of 0.0003 is increased to 0.0009 
or 0.0018 for “R78” or “R508,” respectively). The overall rms and the rms in the interior 
are not much affected, as may be seen from Table 1; however, if relativistic corrections had 
been included in the equation of state, a decrease in the solar radius probably would have 
led to some improvement. 


Solar luminosity effects: Bahcall et al. (2001) tested the effect of 2-a changes in the value 

of L & (namely, ±0.8%) on their solar models, finding only a minor effect on neutrino fluxes 
and negligible effects on the other quantities they considered. For completeness, we made 
the same test with our own models, confirming their results. Figure la shows the effects on 
the sound speed of a solar luminosity 0.8% lower (“Ljoto”) and 0.8% higher ( T/ii 9 /i ) than 
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the most recent value of 3.842 x 10 33 erg s -1 (Bahcall et al. 2001; Frohlich & Lean 1998; 
Croinmelvnck et al. 1996) — note that our reference standard solar model lies closer to the 
latter than to the former, as it uses a slightly higher L© value than the most recent estimate 
(see § 2). To avoid confusion with other curves, the “ Lh tg h ” and “L/ow” curves are shown 
only in the region where they differ the most, namely, r < 0.3 i?©; even in this region, a shift 
of 0.8% in L© produces a fractional change in the sound speed of less than 3 parts in 10 4 , 
dropping to 1 part in 10 4 for r > 0.3 i?©. 


Solar age effects: Figure la demonstrates that the uncertainty in the solar age f© has only 

a very minor effect on the solar sound speed profile — note that the shift illustrated here, 
from t e = 4.6 Gyr to 4.5 Gyr, is much larger than the observationally allowed range of solar 
ages (i.e., 4.55 Gyr < f© < 4.59 Gyr) discussed in § 2. (These ages are defined to include 
the pre-main-sequence; main sequence ages can be obtained by subtracting 0.04 Gyr). The 
maximum allowed shift of 0.02 Gyr relative to the “best” solar age of 4.57 Gyr would yield 
negligibly small effects, namely, an rms of 0.0001 for A c/c and an rms of 0.001 for A p/p. 
with maximum changes less than twice these values. Our results agree both qualitatively 
and quantitatively with the age sensitivity found in the recent work of Morel et al. (1997). 


Low-temperature opacity effects: Uncertainties in the low-temperature molecular opac- 

ities appear to have negligible impact on the structure of a solar model. Using the Sharp 
(1992) molecular opacities (“K S har P ”) below 10 4 K rather than the Alexander k Ferguson 
(1994) molecular opacities ( “/^Alexander” ) led to essentially identical sound speed and density 
profiles -- the “/c S har P ” case is thus not plotted in Figure 2. The rms differences are negligible, 
less than a part in 10 4 for the sound speed and less than a part in 10 3 for the density. 


Interior opacity effects: Recently, Morel et al. (1997) have demonstrated the serious 

impact of opacity changes on the sound speed and density profiles of solar models. They 
compared the 1995 OPAL opacities (Iglesias & Rogers 1996; Rogers et al. 1996) with the 
less-precise 1992 OPAL opacities (Rogers & Iglesias 1992) (albeit with models that neglected 
diffusion), finding that the improved opacities made an improvement of up to 0.005 in the 
sound speed and up to 0.03 in the density (see their models Si and S2). Basu et al. (2000) 
compared a model with the 1995 OPAL opacities and the OPAL equation of state (Rogers 
et al. 1996) to a model with the 1992 OPAL opacities and the cruder Yale equation of 
state (Guenther et al. 1992) with the Debye-Hiickel correction (Bahcall, Bahcall, & Shaviv 
1968) (their models did include diffusion); they likewise found an effect of up to 0.005 in 
the sound speed and up to 0.03 in the density, from the combination of these two changes 
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in the input. We found that an even larger improvement of up to 0.007 in the sound speed 
resulted from changing from the even older Los Alamos (LAOL) opacities (Keady 1985, 
private communication) to the 1995 OPAL opacities, with an mis improvement of 0.004 as 
shown in Table 1. (there is also an improvement of up to 0.04 in the density, with an rms 
improvement of 0.02). However, such a large opacity change overestimates the uncertainty 
in the 1995 OPAL opacities (Rogers k Iglesias 1998). 

As pointed out by Morel et al. (1997), neglecting the opacity changes that result from 
metallicity variations in the Sun would lead to significant errors — - e.g., errors of up to 0.0015 
in the sound speed. In our reference standard solar model, we did our best to account for these 
temporal and spatial variations in the opacity due to these composition changes from diffusion 
and nuclear burning. The metallicity value Z K that we used for metallicity interpolation in 
the OPAL opacity tables was scaled according to the changes in the elements heavier than 
oxygen (resulting from diffusion): Z K — Z h , where Zh = Zo[Ylheavy ^ ie 

extent to which the CNO element abundances differ from the scaled solar values in Z h was 
accounted for by an extra interpolation in nominal “excess carbon and oxygen” abundances 
C ex and O ex using the OPAL opacity tables for extra C and/or O abundances. (Several 
similar methods of doing this were considered; as discussed in § 2, all gave essentially identical 
results — even omitting the CNO-correction entirely had almost no effect.) 

An alternative approximation ( “approx-/c” ) is to set Z K — Z (i.e., to interpolate the 
OPAL opacity tables in the local metallicity Z, but ignore effects of variations in the makeup 
of Z). As shown in Figure 2, this approximation yields results almost identical to those of 
our reference standard solar model; the rms difference was only 0.0002 in A c/c and 0.001 in 
A pj p. Morel et al. (1997) compared two different ways of estimating the value of Z h (both 
being similar but not identical to our “approx-/c” case); they likewise found only very minor 
differences in the sound speed between their two methods, but much larger effects on the 
density (see their models D3 and D12). 

Even when using the most up-to-date OPAL opacities, one can still get significant errors 
if one neglects the effect on opacity of Z-changes (primarily due to diffusion). The simplest 
case is to set Z K = Z 0 ( “const- Z K ”); in effect, such a case uses only the OPAL opacity tables 
relevant to the protosolar metallicity and ignores the effect on the opacity of any subsequent 
changes in the metallicity. Figure 2 demonstrates that this “const- Z*” case yields errors of 
up to ~ 0.0015 in A c/c relative to the more accurate opacity interpolation of the reference 
standard solar model, in agreement with the results of Morel et al. (1997) (compare their 
models D10 and D12); the rms differences were 0.0010 in A c/c and 0.007 in A p/ p. One might 
attempt to fix up this neglect of metallicity variation by interpolation using the mildly CO- 
enhanced OPAL opacity tables, i.e., retaining a constant Z K — Z 0 for opacity interpolation 
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purposes but setting C() fX to the difference between Z K and the true value of Z that results 
from diffusion and nuclear burning, i.e., CO ex — Z — Z K = Z — Z 0 . Figure 2 illustrates 
that this approximation is a slight improvement over the “const- Z K ” case but still not very 
satisfactory: rms differences relative to the reference standard solar model are 0.0006 in A c/c 
and 0.004 in A pj p. 

Interpolation errors can also arise from the finite grid spacing of the opacity tables in 
A', Z, T, and p. In creating the opacity interpolation routines, we tested the A- and Z- 
interpolation, finding that these should result in only minor errors (a fraction of a percent). 

In addition to errors introduced by methods of interpolating in opacity tables, one must 
consider the errors in the actual opacity values contained in the tables. Such opacity errors 
will in general be functions of temperature and density; also, different elements will have 
different errors. Rogers & Iglesias (1998) estimate that there is a 4% uncertainty in the 1995 
OPAL opacities from effects neglected in their calculations. Turcotte et al. (1998) showed 
that differences in the Rosseland mean opacities between 1992 and 1995 OPAL opacities 
do not exceed 7% over the run of temperature and density in the Sun’s interior, and that 
these opacity differences yielded sound speed differences of up to 0.002 in their solar models. 
In addition, errors in the observed relative heavy element abundances in the solar envelope 
will translate into opacity errors, since different elements have somewhat different opacities. 
Rogers & Iglesias (1998) estimate that such abundance uncertainties correspond to opacity 
uncertainties of order 5% at temperatures where ionization effects of the relevant elements 
yield a large contribution to the opacity (e.g., near 2 x 10 6 K for oxygen or neon). Turcotte 
et al. (1998) found that taking into account the opacity effects due to changes in the relative 
abundances of all the individual elements in Z led to opacities that differed by up to 2% 
from the opacities tabulated for the standard scaled-solar metallicity, yielding sound speed 
differences of up to 0.001 and density differences up to 0.005. 

We did not attempt detailed element-by-element variations of the OPAL opacities in 
our models; nor did we test the effect of opacity variations in limited density or temperature 
ranges. Instead, we obtained a rough estimate of the maximum possible effects of uncer- 
tainties in heavy-element opacities by making an overall shift in the metallicity value Z K 
used for interpolation in the OPAL opacity tables. Figure 2 illustrates the cases Z K = 0.9 Z 
(“1ow-k”) and Z K = 1.1 Z (“high-/c”). This shift in the value of Z K corresponds to an average 
shift in the opacities of order 10% over the solar interior (2 — 5% for r < 0.4 R 0 , 10% for 

0.4 Rq < r < 0.7 R q , ~ 15% for 0.7 Rq < r < 0.92 R Q , and ~ 5% for r > 0.92 R e ). As 
may be seen from Figure 2, such an opacity shift of order 10% yields changes of up to 0.003 
in A c/c; the rms differences are 0.0016 in A c/c and 0.014 in A pf p. 

We tested cases where new OPAL opacities had been calculated for mixes in which the 
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abundances of C, N, 0, and Ne were either all increased ( " ^oPAL:GN 93 tc-Nc ,: ) by their quoted 
errors of 15% or decreased by this amount (“/coPAL:GN 934 C-Ne” ), relative to their abundances 
quoted by Grevesse & Sauval (1998). Figure 2 illustrates an effect of up to 0.001 in the 
sound speed from such a change. However, we used a value of ZjX = 0.0245 for all of 
the models in Figure 2, which is not strictly consistent with such large abundance changes: 
since C, N, O, and Ne comprise the major portion of Z , a 15% change in their abundances 
should correspond to a change of ~ 12% in ZjX as well. We discuss such a self-consistent 
composition-plus- .Z’/A'’ variation under ZjX effects below. 


Solar abundance (Z/X) effects: Our reference standard solar model used the obser- 

vational value of Z/X = 0.0245 from Grevesse & Noels (1993), since their mixture was 
the one for which the standard OPAL opacity tables (“/copal:GN 93 ”) were available. The 
close-dashed curve in Figure 3 demonstrates the effects of using the Grevesse & Noels (1993) 
relative metal abundances with the corresponding “kopal:GN 93 ” OPAL opacities, but using a 
13% higher value of Z/X = 0.0277, the older value that had been recommended by Grevesse 
(1984); the maximum sound speed difference relative to the reference standard solar model 
is 0.0030, and the rms differences are 0.0018 in A c/c and 0.017 in A p/p. The dotted curve 
in Figure 3 illustrates a similar case with a 6% lower value of Z/X = 0.023, as recommended 
by the more recent work of Grevesse & Sauval (1998); the maximum sound speed difference 
is 0.0016, with rms differences of 0.0010 in A c/c and 0.009 in A p/p. (Note that most of the 
above effect comes from the different opacity that results from the changed solar Z value, 
as may be seen by comparing with the “high-Ac” and “Iow-k,” curves in Fig. 2.) However, 
these comparisons are not strictly self-consistent, since it is the changes in the individual 
elemental abundances of the metals that adds up to yield the changed ZjX ratio. Using the 
old abundance pattern of Grevesse (1984) and newly-computed OPAL opacities appropriate 
to it (“/coPAL:Gr 84 ”) leads to the wide-dashed Z/X = 0.0277 curve in Figure 3, with a slightly 
smaller maximum sound speed difference (of 0.0025) and rms differences (of 0.0014 in A c/c 
and 0.011 in A p/p). Similarly, using the Grevesse & Sauval (1998) abundance pattern and 
appropriate OPAL opacities (“«opal:GS98”) leads to the thin solid Z/X = 0.023 curve in 
Figure 3, reducing the maximum sound speed difference to 0.0012 and the rms differences 
to 0.0006 in A c/c and 0.004 in A p/p. 

Strictly, the uncertainty resulting from observational solar abundance errors can be 
estimated by varying the solar abundance values of Grevesse & Sauval (1998) within their 
quoted uncertainties, obtaining OPAL opacities with these revised compositions, calculating 
the resulting Z/X values, and running solar models with these self-consistent sets of input 
values. We have done this for two cases. Rather than performing large numbers of random 
variations of the abundances, we tested a case which should give something close to the 
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maximum effect. The elements C, N, O, and Ne not only comprise the major part of the 
metallicity but also have relatively large errors of ~ 15%, and unlike other elements with 
large errors one cannot get a “better” value by using the meteoritic abundance instead. We 
therefore considered cases where C, N, O, and Ne were either all increased by 15% ( Z/X = 

0. 0257, “/C0PAL.GS98TC-Ne”) or all decreased by 15% (Z/X = 0.0203, “/copAL:GS98|C-Ne” ) — 

1. e., these self-consistent abundance variations correspond to 12% variations in Z/X . As 
illustrated in Figure 3, these cases lead to variations in the solar sound speed of up to 0.003 
relative to the Z/X = 0.023, “«opal:GS 98 ” case, with rrns differences of 0.0017 in A c/c and 
0.012 in A p/ p. 

Morel et al. (1997) considered the effect of a 6% increase in Z/X , finding a maximum 
difference of 0.0007 in their sound speed and 0.003 in their density; this would imply a 
significantly lower sensitivity to Z/X than we found. This is probably due to the fact 
that the models in which they tested Z/X variations did not consider the effect on the 
opacities of the temporal and spatial variations in the heavy element abundances that arise 
from diffusion, but merely used opacities appropriate to a constant metallicity equal to the 
protosolar value ( Z K — Z 0 , as in our “const- Z K ” case discussed above). 


Nuclear rate effects: Figure 4 demonstrates that the uncertainty in the ] H(p, ive + ) 2 H 

nuclear burning rate (the pp rate) has a significant impact on the solar sound speed. Our 
reference standard solar model used the recommended nuclear rates from the NACRE com- 
pilation (Angulo et al. 1999). These authors also supply “high” and “low” cases to indicate 
the allowed uncertainty range of each nuclear rate; in the case of the pp reaction, the high 
case is 8% above the recommended rate and the low case is 3% below it. We have tested 
the effects of nuclear rate uncertainties by computing variant standard solar models using 
high and low NACRE rate values. Figure 4 demonstrates that a high pp rate is preferable, 
if all other parameters are kept constant: the high pp rate gives good agreement with the 
helioseismic reference profiles, except in the Sun’s central regions where the helioseismic ob- 
servations are the poorest. Our models indicate that a change of 5% in the pp rate yields 
changes of up to 0.003 in the sound speed, (0.0014 in the regions accurately probed by he- 
lioseismology, outside the core); the rms changes in such a case would be 0.0009 in A c/c and 
0.018 in A pj p. Antia h Chitre (1999) also tested the effects of changes in the pp rate on the 
helioseismic profiles, concluding that a relatively high pp rate is preferred, consistent with 
our results discussed above. 

Figure 4 also demonstrates that the uncertainty of ±20% in the 3 He(er, 7 ) 7 Be reaction 
leads to only a minor effect: a maximum change of 0.001 (or 0.0003 outside the core), with an 
rms relative change of 0.0002 in A c/c and 0.004 in A p/p. Basu et al. (2000) considered the 
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effect of setting 3 He (a, 7) 7 Be rate to zero; they found large effects from such an unphysically 
extreme change. Setting the rate to zero is equivalent to a 100% change, 5 times as large as 
the 20% considered by us; thus it is consistent that their published effect is about 5 times 
as large as ours. 

We also tested the effects of the ±6% uncertainty in the 3 He( 3 He, 2p) 4 He reaction and of 
the ±30% uncertainty in the 14 N(p,7) 15 0 reaction (which determines the CNO-cvcle rate). 
Such changes in these rates led to negligible effects on the sound speed and density profiles; 
we have not plotted these profiles in Figure 4, since they would be essentially superimposed 
on that of the reference standard solar model. 

It is not surprising that the uncertainty in the pp rate has the largest effect on the sound 
speed and density profiles, since it is the basic rate that determines the overall p-p chain 
burning rate. 

We also computed a model using the previous standard set of nuclear rates, namely, the 
Kellogg nuclear rate compilation of Caughlan k Fowler (1988). The resulting sound speed 
and density profiles are shown in Figure 4. With us still in Kellogg, carrying out this work in 
an office directly below his long-time office, it is especially gratifying for us to see that Willy 
Fowler’s last published pp rate yields such good agreement with the current helioseismic 
reference profiles (the largest differences being near the center, where the observations are 
least accurate). 


Diffusion effects: There are uncertainties in the diffusion coefficients; Proffitt (1994) es- 

timates a 15% uncertainty in the diffusion constant of helium relative to hydrogen, and a 
~ 50% uncertainty in the diffusion constant for oxygen relative to hydrogen. Figure 5 shows 
that an increase or decrease of 20% in the helium diffusion constants has only a modest effect: 
a maximum change of 0.001 in the sound speed, with rms changes of 0.0008 in Ac / c and 
0.007 in A p/p. The effect of increasing or decreasing the heavy element diffusion constants 
by 40% has an even smaller effect: a maximum change of 0.0006 in the sound speed, with 
rms changes of 0.0004 in A c/c and 0.004 in A p/p. 


3.2. Solar Convective Envelope Depth 

One of the key results of helioseismic observations is a highly precise value for the 
position Ra, of the base of the solar surface convective region: Basu k Antia (1997) report 
a value of R ce = 0.713 ± 0.001 Rq. Our reference standard solar model is in agreement with 
this value, having R ce = 0.7134 Rq (see Table 1), independent of whether we used the OPAL 
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or the MHD equation of state at log p < -2. The low and high Z/X values are marginally 
inconsistent with the observed value, lying at R ce — 0.7181 and 0.7096 Rq, respectively. 
The uncertainty of ±0.05 Gyr does not make a significant difference in R ce (see Table 1). 
As far as the uncertainties in the nuclear reaction rates are concerned, only the pp rate 
has a significant effect on R ce , of ±0.002 Rq. Uncertainties in molecular opacities do not 
have a significant effect on R ce , but changing from the old 1985 LAOL opacities to the 1995 
OPAL opacities does yield a large improvement (of 0.006 /?©); the remaining uncertainties 
in the 1995 OPAL opacities might thus be expected to have a small but possibly significant 
effect on R ce . The uncertainty in the diffusion constants for helium does have a significant 
influence (±0.003 R Q ), but uncertainties in the diffusion constants for the heavy elements 
do not (effects < 0.001 Rq). 

Note that the cases favored by the sound speed profiles at the one- to two-sigma signif- 
icance level (high opacities, high Z/X, or high pp rate) are disfavored by the observed R ce 
value at about the same significance level. 


3.3. Solar Helium Abundance 

Another key result of lielioseismic observations is a precise value of solar envelope helium 
mass fraction Y e (this value is lower than the Sun’s initial helium abundance, due to diffu- 
sion). Inferring the solar helium abundance requires use of a (theoretical) equation of state 
as well as the lielioseismic frequency observations that probe the solar convective region, 
particularly the He II ionization zone at 0.975 Rq < r < 0-985 Rq (Richard et al. 1998). In 
recent work, there has been excellent agreement on the value of Y e when the OPAL equation 
of state was used: Y e = 0.249 ± 0.001 (Basu & Antia 1995), Y e = 0.248 ± 0.006 (Kosovichev 
1997), Y e = 0.248 ± 0.001 (Basu 1998), and Y e = 0.248 ± 0.002 (Richard et al. 1998). A 
larger scatter has resulted when the MHD equation of state was used: Y f . - 0.242 ± 0.003 
(Perez Hernandez & Christensen- Dalsgaard 1994), Y e = 0.246 ±0.001 (Basu & Antia 1995), 
i; = 0.232 ± 0.006 (Kosovichev 1997), Y e « 0.252 (Basu 1998), and Y e ^ 0.242 (Richard 
et al. 1998). The OPAL equation of state is expected to be more accurate than the MHD 
equation of state over the bulk of the convective envelope, since it was designed for use at 
log p < -2 (which occurs in the Sun at r > 0.942 R Q ); however, the He II ionization zone 
occurs in a region the MHD equation of state was specifically designed for, and w'here it may 
actually be more accurate than the OPAL equation of state (see, e.g., Richard et al. 1998). 
In summary, the helioseismic observations appear to constrain the solar envelope helium 
abundance to the range 0.24 < Y e < 0.25. 

Our theoretical reference standard solar model is in excellent agreement with this, having 
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Y e = 0.2424 independent of whether we used the OPAL or the MHD equation ot state at 
log p < -2. Low and high Z/X values yield Y e values of 0.2344 and 0.2517, respectively, 
on the borderline of acceptability. Uncertainties in the solar age, the nuclear rates, and 
opacities have only a negligible effect on Y e \ uncertainties in the diffusion constants have a 
minor effect, but still yield acceptable values (see Table 1). 


3.4. Solar Lithium Abundance 

The present observed solar surface lithium abundance is loge('Li) = 1.10 ± 0.10 as 
compared to the initial value of log e( 7 Li) = 3.31 ± 0.04 obtained from meteorites (Grevesse 
& Sauval 1998), where log e( 7 Li) = log (Ay/An) + 12 for number densities Ay and Ah of 
lithium and hydrogen, respectively. The solar surface lithium depletion factor / Li , relative to 
its initial value, is thus observed to be /y = 160 ± 40. Solar surface lithium can be depleted 
due to three causes: (1) lithium burning during the pre-main-sequence evolution, when the 
surface convection still reaches deeply into the interior: (2) rotationallv induced mixing on 
the main sequence, which transports lithium down from the convective envelope to regions 
hot enough for lithium burning; (3) mass loss on the main sequence, which can cause the 
convective envelope to move inwards and engulf lithium-depleted regions. In this paper, we 
only consider the first of these, namely, the pre-main-sequence lithium destruction; rotational 
mixing is beyond the scope of this paper, and main sequence mass loss is discussed in the 
companion “Our Sun V” paper (Saekmann & Boothroyd 2001). Our reference standard solar 
model had a pre-main-sequence lithium depletion factor / Li = 24, as shown in Table 1. 

Pre-main-sequence lithium depletion depends quite sensitively on the structure of the 
solar models during that stage of evolution. We give our results in terms of the solar lithium 
depletion factor /y, which is shown in Table 1. The zoning did not affect the lithium 
depletion factor significantly; nor did the uncertainties in the solar age, luminosity, and 
radius. Changes in the equation of state in the outermost regions can have a small effect 
(< 30%) on the lithium depletion factor: the “OPALeos-midT” and “OPALeos-lowT” cases 
had lithium depletion factors of f u = 17 and 19, respectively, as opposed to f u = 24 for the 
standard “OPALeos-hiT” case. 

The pre-main-sequence lithium depletion is extremely sensitive to both low-temperature 
and high-temperature opacities; use of the Sharp (1992) molecular opacities instead of the 
Alexander & Ferguson (1994) ones halved the lithium depletion factor (/y = 10), as did our 
“1ow-k” test case (f u = 10), while our “high-*” test case nearly tripled it (/y = 71). There 
is also a relatively large sensitivity to the uncertainty in the observed solar abundances. A 
low Z/X ratio of 0.0203 (with U K 0 PAL:GS 984 C-Ne”) yields only two-thirds as much lithium 
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depletion (fu — 15), while a high Z/X ratio of 0.0257 (with “KoPAL:GS>iKte-N<‘ ) yields hall 
again as much lithium depletion (fu = 35). 

Uncertainties in the diffusion constants can affect the lithium depletion factor signifi- 
cantly, primarily due to the effect of the different initial composition; our diffusion test cases 
have depletion factors ranging from fu — 18 to 33. 

Except for the 7 Li + p rate, uncertainties in the nuclear rates have almost no effect on 
the extent of lithium depletion. For the 7 Li ± p rate, the ±14% uncertainty quoted by the 
NACRE compilation (Angulo et al. 1999) corresponds to an uncertainty of about ±50% in 
the depletion factor (i.e., a range in the from fu = 16 to 38). 


Solar beryllium abundance: The observed solar beryllium abundance is log e( 9 Be) = 

1.40 ± 0.09, consistent with no depletion relative to the meteoritic value of log e( 9 Be) = 
1.42 ± 0.04. These values imply that solar beryllium cannot have been depleted by more 
than a factor of 2 (3-<7 upper limit). Our solar models all had negligible amounts of beryllium 
depletion, of order 1%. 


3.5. Solar neutrinos 

We will not devote much space to the predicted solar neutrino values, since it has 
long been concluded that matching the observed neutrino capture rates requires not revised 
astrophysics but new neutrino physics, e.g., Mikheyev-Smirnov-Wolfenstein (MSW) neutrino 
oscillation effects (see, e.g., Bahcall, Basu, & Pinsonneault 1998a; Suzuki 1998; Basu et 
al. 2000; Bahcall et al. 2001). In Table 1 we present the theoretically predicted neutrino 
capture rates for the 37 C1 and 71 Ga experiments, and the predicted flux of 8 B neutrinos. 
As is normally obtained, our theoretical predicted neutrino rates are much in excess of the 
observed values, i.e., 6.4 to 8.9 SNU is predicted for the :,7 C1 experiment, as compared to the 
observed value of 2.56 ± 0.23 SNU (Davis 1994; Cleveland et al. 1998; Lande 2001), and 127 
to 141 SNU for the 71 Ga experiments, as compared to the observed value of 74.5 ± 5.7 SNU 
(combined value from SAGE and GALLEX±GNO: Hampel et al. 1999; Abdurashitov et 
al. 1999; Altmann et al. 1999; Gavrin 2001; Ferrari 2001). Likewise, the models predict 8 B 
fluxes of 4.4 to 6.3 x 10 6 cm -2 s _1 , as compared to the value of (2.32 ± 0.08) x 10 6 cm -2 s -1 
measured by Super-Kamiokande (Fukuda et al. 2001). 
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Fig. 1. — Changing the zoning, the outer-envelope equation of state, or the solar radius, 
luminosity, or age: the effects on (a) the adiabatic sound speed c, and (b) the adiabatic in- 
dex r, . The fine-zoned reference standard solar model (solid line — errorbars give statistical 
error in the inferred helioseismic profile) switches from the OPAL to the MHD equation of 
state as the density drops below log p = —2 (i.e., r « 0.94 R & )\ the coarse zoned ( dotted 
line) and very coarse zoned ( dot-dashed line) models are the same except lor the zoning. 1 he 
fine-zoned “OPALeos-midT” ( short-dashed line) and “OPALeos-lowT ’ ( dashed line) models 
have the equation-of-state switchover at log T — 4.0 and log T = 3.75, respectively. The 
wide-dashed line shows the effect of using a reduced solar age (actually, ~ 10-er below the 
lower limit of 4.55 Gyr). For clarity, luminosity and radius effects are shown only in regions 
where they have a visible effect. For r < 0.3 Rq , thin solid lines show the effects of using 
the maximum (“ Lhigh ”) and minimum (“L( olu ”) values of Lq. For r > 0.71 /?©, thin solid 
lines show the effects of using R & = 695.78 Mm (“R78”) or R & = 695.508 Mm (“R508”), 
while thin dashed lines show the same radius effects for the “OPALeos-midT” case. 
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Fig. 2.— Effect on the sound speed of opacity uncertainties. The dot-dashed lines (“low- 
K r and “high-/c”) show the effect of a ~ 10% overall change in the opacities relat ive to the 
reference standard solar model ( thick solid line). Ignoring opacity effects due to relative CNO 
abundance changes from CNO-burning has little effect ( dashed line), but ignoring opacity 
effects due to metallicity variations from diffusion leads to large errors {short- dashed and 
wide-dashed lines). Opacity uncertainties from uncertainties in relative metal abundances 
can be significant, even when the solar surface ZjX ratio is unchanged (thin solid and dotted, 
lines). 
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Fig. 3. — Effect on the sound speed of uncertainties in the observed solar surface composition. 
Relative to the reference standard solar model at Z/X = 0.0245 ( thick solid line), switching 
to the more recent Z/X = 0.023 value of Grevesse & Sauval (1998) without including the 
corresponding changes in relative metal abundances ( dotted line) has a larger effect than the 
full case with self-consistent abundances and opacities (‘‘/copal:GS98 ; ’ thin solid line); the 
same is true when considering the older Z/X — 0.0277 cases of Grevesse (1984) (dashed line 
vs. self-consistent “/coPAL:Gr 84 ” wide-dashed line). Nonetheless, a self-consistent estimate of 
the uncertainties in the Grevesse & Sauval (1998) composition ( dot-dashed lines) shows a 
large effect (compare to the thin solid line). 
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Fig. 4. — Effect on the sound speed of uncertainties in nuclear rates. The ~ 5% uncertainty 
in the basic pp rate (dot-dashed lines) has a much larger effect than the 20% uncertainty in 
the 3 He + a rate ( dotted lines)-, other rate uncertainties have negligible effects, and are not 
plotted. The difference between the reference standard solar model using the NACRE rates 
( solid line) and the model using the Caughlan & Fowler (1988) rates (“CF88:” dashed line) 
is largely due to the higher pp rate adopted by the latter authors. 



Fig. 5. — Effect on the sound speed of uncertainties in diffusion constants. A 20% variation 
in the diffusion constant for helium ( short-dashed and long-dashed lines) has a larger effect 
than a 40% variation in the heavy element diffusion constants ( dotted and dot-dashed lines) 
— the reference standard solar model is shown by the solid line. 





Table 1. Characteristics of Our Solar Models' 
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g For these cases, the relative rms values compare the variant models with the case relative to which the parameter variation was performed. 
h Most recent solar luminosity value of L& = Lb est = 3.842 x 10 3 ^ erg s _1 . 

1 OPAL opacities not interpolated in “excess** C and 0. i.e.. opacities computed as if the C’NO clement abundance profiles in the Sun were always in the 
same proportion to the Fe profile. 




